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CHAPTER I 


INTRODUCTION 


1.1 General Statement of the Problem 
When a nuclear weapon explodes, a high pressure wave develops 

and moves away, from the point of the explosion, traveling over the ground 
surface and through the ground. A _ shallow-buried arch subjected to this 
pressure responds in various ways, depending on the orientation of the 
structure relative to the direction of travel of the blast wave, the parame- 
ters of the blast loading, and the properties of the soil surrounding the 
arch. If the arch is oriented with its longitudinal axis perpendicular to 
the direction of travel of the blast wave, the windward side of the arch 
is subjected to the pressure a few milliseconds before the leeward side. 
This unsymmetrical sequence of loading produces flexural stresses and de- 
flection in the arch. The deflection of the arch into the surrounding soil 
gives rise to an unsymmetrical passive earth pressure loading. 

; An understanding of how these loadings occur and their effects on 


the arch is necessary for the proper consideration of the design of the 


arch rib. 


Vee Purpose of this Thesis 


When the arch rib is subjected to the above described loads it 
responds in a configuration which is a function of an infinite number of 
shapes or modes of vibration, since the arch rib has distributed mass and 
elasticity. An infinitely large number of coordinates would be required 
to specify the position of the deflected arch. The participation of each 


mode depends on the spatial distribution of the loads along the arch rib 





2 
as well as the time variation of the dynamic load in relation to the natural 
frequencies of vibration of each mode. 

The distortion of the arch rib is a function of the stiffness of 
the rib, the passive earth pressure loading resulting from the stiffness of 
the soil surrounding the arch, the distribution of the soil and dynamic 
loads, the time variation of the dynamic load, and the support condition of 
the abutments. It is convenient to consider the earth pressure loading due 
to the stiffness of the soil as an integral part of the resistance of the 
arch. For flexible arches, the stiffness of the soil is greatly in excess 
of the stiffness of the arch rib, consequently the properties of the arch 
rib may be neglected. Henceforth, unless specifically designated otherwise, 
the resistance of the arch rib means the resistance given to the arch by 
the soil. 

Design of structures for dynamic loading is usually based on the 
principle of conservation of energy. When the dynamic load acts on the 
structural element, the element deflects. Work is done on the structure by 
the load at any point as that point deflects. Due to the speed of the load- 
ing, kinetic energy of motion is given to the mass of the element, and, as 
the element deflects, strain energy is stored. The strain energy may involve 
both elastic and plastic action. A general statement of the principle of 
. conservation of energy is: 

work done = kinetic energy + strain energy 

To evaluate this equation the displacement of any point on the 
structure under the load must be known. The strain energy is the sum of 
all the strain energy inthe element, and may be caused by compressive, 
flexural, torsional and shear stresses in the element. The kinetic energy 


is the energy of translation or rotation of all the masses in the element. 





-_ 
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The evaluation of these quantities would be exceedingly difficult 
if no more than the arch rib were involved. Evaluation of the energy 
stored in the soil is a hopeless task, impossible for practical purposes. 
Fortunately, the principal modes of response ofr the arch rib are a syn- 
metrical mode corresponding to the compression stresses and an antisymmetrical 
mode corresponding to the flexural stresses. This allows a considerable 
simplification to be made in the treatment of the problen. 

A single-degree-of-freedom system consists of a lumped mass 
restrained by a weightless spring. When this system is subjected to a load, 
only one coordinate is required to fully define the position of the mass. 

By assuming that the antisymmetrical mode of response of the arch rib is 
defined by the radial deflection of a point at the quarterpoint of the arch, 
which is the point of maximum deflection of the arch, the entire deflected 
shape of the arch may be defined with one coordinate, and the arch may be 
considered a single-degree-of-freedom system. In essence this neglects any 
tangential movements of the arch rib associated with the radial deflection. 

After representing the action of the arch rib as that of a single- 
degree-of-freedom system, dynamic equivalence relations should be established 
between the actual and idealized masses, spring stiffnesses, and loads acting 
on the two systems. However, for the purposes of this study, this refinement 
is omitted. Instead, the following assumptions are made: 

(1) The mass of the single-degree-of-freedom system is equal to 
the unit mass of the arch, which is obtained by uniformly distributing along 
the Bch rib the weight of the soil which creates the passive earth pressure. 

(2) The resistance function of the single-degree-of- freedom 


system is equal to the unit resistance given the arch by the earth pressure. 





(3) The forcing function on the single-degree-of-freedom system 
is equal to the unit pressure of the forcing function acting on the arch rib. 
Once the mass, resistance function, and forcing function acting 
on the single-degree-of-freedom system have been defined, the response of 
the system may be determined by the well-known differential equation of 
motion of such‘a systen: 


mx + kx = p(t) 


where x = the second derivative of the displacement with respect 
to time 
kx = the resistance of the system as a function of the 
displacement 
m = mass Of the system 


p(t) the externally applied force as a function of time 


The main purpose of this thesis is to investigate the response of 
various arch ribs under a limited range of variables when acted upon by blast 
pressure. In order to do this, it is necessary to relate the parameters of 
the single-degree-of-freedom system to those of the actual structure under 


investigation. 


1.4 Method of Approach 


The blast pressure acting on the surface of the ground causes 
vertical and horizontal pressures in the soil. These pressures are a 
function of the maximum surface pressure, and are time-dependent upon a 
pressure wave which passes through the soil. Pressure-time relations of 
these pressures are formulated. 

These pressures are assumed to be principal stresses, the resultant 
of which acts radially on the arch. The tangential component is omitted. 


The total radial load which results from these pressures is designated as 





a load profile PA? t), which is then represented as a distributed load by 


the series ) p(t) sin == 





2 XK 
where p(t) a [- p_.(, t) gin = Y Ra oy) 
O 


@ = central angle of the arch 


R = arch radius 


An analogous, uniformly-distributed, antisymmetrical, time- 


dependent load, p alt): is then determined. This load has the same general 





effect on the entire haunch of the arch as the 1oea ) P, (t) sin = zy. 
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The pressure-time diagram for p(t) and the numerical relation between P( t) 
and the dynamic load) p(t) sin =< 2 are established by equating their 


effects at various times as though they both acted statically. The analogous 
antisymmetrical load p(t) defines the forcing function acting on the arch 
and also on the single-degree-of-freedom systen. 

The unit, horizontal, passive earth pressure acting on the arch 
is determined by Rankine's theory for passive earth pressure. The unit 
vertical pressure acting on the arch is found by statics. These two pressures 
are assumed to be principal stresses, and their radial resultant is assumed 
to resist the movement of the arch. The tangential component is omitted. 
The total radial soil load is designated as a load profile P(?), and is 


represented as a distributed load by the series ) P,, sin ~—— 


@ 
where p= =e p,.(?) sin 222 Rag 
° 





n OR ) 


An analogous, uniformly-distributed, antisymmetrical load which 


has the same general effect on the entire haunch of the arch as the load 








@ 
assumed to reach its maximum value when a point located at @/4 on the arch 


) Pa sin Bree is determined by statics. The passive earth pressure is 
rib deflects radially an amount equal to 0.001R. In addition, the resistance 
diagram of this pressure is assumed to be elasto-plastic. The analogous 
antisymmetrical load defines the unit resistance of the arch rib and also 
the resistance function for the single-degree-of-freedom system. 
The weight of the soil that causes the passive earth pressure is 

determined by Rankine's theory and by statics. This weight is assumed to 
be uniformly distributed along the arch rib and defines the unit mass of the 
arch rib and also the mass of the single-degree-of-freedom system. 

The response of the equivalent single-degree-of-freedom system is 
determined by a method developed by Dr. N. M. Newmark for blast resistant 


design, Reference (8). 





CHAPTER II 


The variables of the problem may be divided into three general 
groups corresponding to the blast loading, the earth pressure loading, and 


the configuration of the arch rib. 


2.1 Blast Pressure Loading 
The blast forces acting on the arch rib are conceived to act as a 


pressure wave passing through the soil, Figure 2.1. The direction and 
velocity of the propagation of this wave in the soil is assumed to be related 
to the shock velocity of the surface pressure wave, U, and to the seismic 
velocity of the soil surrounding the arch, C5: The magnitude of C, is 
assumed to be less than U. 

Within the stressed-soil zone, the vertical pressure p(t) is 
assumed to be equal to the surface overpressure P(t), but to have a dif- 
ferent pressure-time curve to allow for the retardation effects on the shock 
Peenitted through the soil. The horizontal pressure p(t) is assumed to 
be proportional to the vertical pressure. 

The time of arrival a of the pressure wave at any point in the 
soil, or on the arch, is determined by the propagation velocity of the 
wave. The time of rise t. of the vertical pressure to its maximum value Pa 


is assumed to be a function of the vertical transit time t, of the pressure 


t 
wave. The rate of decay of the earth pressure wave is assumed to be the 
same as that for the surface pressure. These relations are shown 
qualitatively in Figure 2.2. 

Study of Figure 2.1 indicates that as the shock velocity of the 


blast wave increases relative to the seismic velocity of the soil, the 








unsymmetrical loading decreases, since the pressure wave approaches the 
arch more nearly vertically. For the usual range of these variables the 
earth pressure wave qualitatively approaches the arch as shown in Figure e.l. 

It should be noted that this pressure wave is considered to have 
been induced solely by the surface pressure. The direct ground shock wave 
associated with’ a nuclear weapon explosion is not considered in this thesis. 
Reference (1) states that for pressures of less than 300 psi and for depths 
of less than fifty feet the direct ground shock may be ignored. 

The arch is assumed to be acted upon by the radial resultant of 
the vertical and horizontal pressures. If the ratio of the horizontal 
pressure to the vertical pressure is as great as 1:1, the radial pressure 
is equal to the vertical pressure. Moreover, the total radial force in- 
creases as the ratio of the horizontal to the vertical pressure increases. 

The effect of the rise time may be seen in Figure 2.2. The 
pressure wave arrives at some point ?) on the arch, Figure 2.1, and the 
pressure increases to its maximum value in some time ti: At some adjacent 
point ?> the pressure is at a different level, depending on the relative 
arrival and rise times. Reference (2) recommends that the rise time be 
taken equal to one-half the vertical transit time. 

The following range of variables will be studied in this paper: 

Ex = 100 and 200 psi with U = 3000 and 4000 fps 
C = 1000 and 200 fps 


S 


t. = t, and 1/2 t, 


ul 


p/p, = 0-25, 0.50, 0.75 and 1.00 


c.c. Passive Earth Pressure Loading 
The unsymmetrical blast load on the arch causes the arch to deflect 


into the soil on the leeward side. This deflection brings into action the 





passive earth pressure of the soil as well as the inertia of the mass of 
soil associated with this pressure, since the blast load acts dynamically. 
This pressure contributes materially to the resistance of the arch and mst 
be considered. 

While precise evaluation of the passive earth pressure under static 
conditions is complex, evaluation under dynamic conditions is equally, or 
even more, complex. The passive earth pressure can vary in the length of 
the arch due either to the properties of the soil or to the construction 
technique used in placing the backfill. Rankine'’s theory has been chosen as 
the simplest method for evaluating this passive earth pressure. 

In the application of Rankine's formulas, Reference (3), the unit 
weight of the soil is assumed to be 120 pounds per cubic foot and the angle 


of internal friction is assumed to be thirty degrees. 


2.5 Structures Considered 

The blast loading on arches with central angles of 180 and 120 
degrees and radii of 150 and 300 inches is investigated to determine the 
effect of the soil variables, arch geometry, and relative rise time on the 
arch loading, and to determine a conventionalized dynamic loading function. 
Subsequently investigated is the response of arches with central angles 


of 180 and 120 degrees and radii of 200 to 600 inches to this conventionalized 


dynamic loading function. 
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Figure 2.2 Vertical Pressure-Time Diagram for ? 





CHAPTER IIT 


GENERAL RELATIONS COMMON TO EACH LOADING FUNCTION 


To predict the effect of the loads on the arch rib, the loads must 
be defined in a convenient manner. 

Certain general relations may be derived for defining either the 
blast load or the earth pressure load. These general relations concern the 
transformation of stress with direction, the method used to determine the 


value of the coefficients p(t) and p, for the distributed radial load 








g @ 
load on the arch, and the basis for assuming analogous loads to represent 


functions ) p(t) sin PLY ana >. p,, sin na which represent the actual 


these distributed loads. 


For purposes of dynamic analysis, the arch is idealized to a 


nx na 9? 
~“e 


is replaced by a uniformly-distributed, antisymmetrical, time-dependent load 


single-degree-of-freedom system. The time-dependent load) P, (t) sin 





p(t) » and the distributed earth pressure load) P. sin = Y is replaced 
by a uniformly-distributed antisymmetrical load Pee These antisymmetrical 
loads act outwardly on one half of the arch, and they act inwardly on the 
other half. 

The direction of the analogous earth load is opposite to that of 
the analogous blast load. These antisymmetrical loads define the effects 
of the actual loads on the arch haunch. 

The values of these analogous loads are found by the application 


of the principles of statics. The value of the analogous load for the soil 


pressure may be found by equating the moment produced by the distributed load 





K 
» Pp, sin = Y at a point 0/4 to those of the analogous antisymmetrical load 





“om 
2a at the same point. The replacement of the dynamic load > Pat*) SS - 2 


with an analogous load requires that a basic assumption be made. It is 





assumed that if two generally similar static loads produce equal static 
effects on a structure at a specified point, then the same two loads applied 
dynamically with similar pressure-time diagrams of equal duration would 
cause equal dynamic effects at the same point. Of course the dynamic 
effects of the loads are not defined. The point considered is 6/4. 

To determine the analogous loads, the moments caused by the anti- 


symmetrical loads p(t) and p, are equated with the moments caused by the 
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addition, the effect of a uniform compressive load on the central third of 


distributed loads ) p(t) sin 2 a! and » pe sin 2 w respectively. In 


the arch coupled with a uniform outward load on the outer thirds of the arch 
is determined. This load provides an analogous load for study of the crown 


of the arch. 


5.1 Stress Transformation with Direction 

Let it be assumed that a horizontal unit stress Py and a vertical 
unit stress Pp, are acting on a segment Rd? of a semicircular arch rib located 
at the angle ? measured from the right abutment. Assuming that P, and Py 
are principal stresses, the radial pressure p.(?) at @ may be found by the 
well-known method for determining the relations between stresses at a point 
on different planes passing through the point, Reference (4). 

Rd? may be assumed to be a plane surface for an infinitesimal 
segment. The vertical projection of the plane Rd? is Rd® cos 9; the hori- 
zontal projection is Rd? sin 9. Equating the sum of the radial forces to 


zero yields the following: 
p_(9) Rdg = Rd cos“@ + p_ Rd? sin’? 
r Ph Vv 


Using the trigonometric identities 


cos“? = {1 + cos 29) 


sin? = 3{1 - cos 29) 
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and cancelling Rd? from both sides of Equation (a), 
) } 
p,.(9) = 3(1 + cos 29) Pat ={1 - cos 29) PY 

= 2p, +P) + oP, - P,) cos 2p (1) 

The ratios of p/P, to be studied are 0.25, 0.50, 0.75 and 1.0. 

In addition, these pressures are time-dependent, so p,.(?) must be a time- 
dependent pressure and is therefore expressed as P(t) when the pressures 


considered arise from the blast pressure acting on the surface of the ground. 


Let p(t) = Cpi(t). Then Equation (1) becomes: 


p(%t) = C p(t) + p(t) (1 - ¢) sino (2) 


This equation is evaluated for various values of C as follows: 


For C = 0.25, p_(#,t) = @.25 + 0.75 sing) p(t) (3) 
For C = 0.50, P(t) =(.50 + 0.50 sin“9) p(t) (4) 
For C = 0.75, p(%,t) =(0.75 + 0.25 sin9) p(t) (5) 
For C = 1.00, p.(9,t) = p(t) (6) 


Equations (3) through (6) express the time-dependent unit radial 
pressure acting on the arch rib at any point 9 for the various assumed 
ratios of the horizontal pressure p(t) to the vertical pressure p(t). 


Evaluation of these equations is set forth in Table I. 


3.2 Distributed Load Functions 

Any radial load profile p (9); such as that shown in Figure 3.1, 
when acting upon a segment of a circular arch rib with a central angle @, 
may be represented as a distributed load by a trigonometric series in the 


following manner: 
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~ a en P . nx? 
p(?) = p, sin a” + Pp Sin == + + p. singe tigeae 





(7) 
The coefficient pe is obtained by multiplying both sides of 


w and integrating over the length of the arch rib. 





Equation (7) by sin ne 


6g 
e np s mp nap 
, P,(9) sin = oF Rd @ - | Py sin 7 sin = Rd? 
oO re) 
¥ emp np 
+ p,, sin —— sin — Rd? 
S 6 6 
fe) 
6 
+ so0- +f DP, siné a Rd? + -- 
o (8) 
2 mrp nx 
The value of the integral | sin —> sin = Rd? is zero provided 
oO 


that m and n are integers and m does not equal n. Equation (8) then yields: 








8 nn , sin 2. 6 
i p_.(?) sin — > Rd? = R py | $- ae 
fe) 5 P 
= 4 P. R6@ 
a : np 
. P= [ p.(?) sin —> Rdp (9) 
fe) 


If p.(%) can be expressed analytically, the value of p_ may be 
found by integration. If p_.(?) is an irregular load profile for which 


there is no convenient analytical expression, the integral may be evaluated 


by numerical methods since the integral f p.(?) sin ae Rad? represents the 


r 
area under the load profile multiplied by sin = 


The sign convention adopted throughout this thesis is that a 


positive load acts inwardly and a negative load acts outwardly. The 
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harmonic components of loading P., sin -- shown in Figure 3.1 are in 


qualitative agreement with the original load p (9) 3; therefore the sign of 
the coefficient Pr is definitely negative. The sign of the coefficient isk 
has been taken as negative, although no positive statement can be made that 
this is the case for the original load p_.() in Figure 3.1. 

Included in Figure 3.1 are two cases of uniform radial loading. 
It can be shown that the effects of) P., sin = on the crown and haunch 
of the arch may be represented by such uniformly distributed loads. 

If the load profile p(?) is time-dependent, it is expressed as 
P(?, t) . The procedure described above may be used to determine the co- 
efficient p(t) ,» with the qualification that p(t) must be evaluated at a 
specified time. 

Reference (5) eenene a full development of the methods of 


harmonic analysis as applied to structural engineering problems. 


5.5. Basis for the Determination of Analogous loads 
As stated at the beginning of this chapter, the analogous loads 


p,(t) and p, are derived by equating their moments to the sum of the moments 


of the distributed loads a p(t) sin a and ) 2 sin -S , respectively. 


The time dependency of p(t) is not considered in this section since the 


load ) P(t) sin = is assumed to be acting statically at the instant 


the analogous load p(t) is evaluated. 
Moment equations for all load profiles in Figure 3.1 are found 
in this section after general equations for the right reaction of any two- 


hinged semicircular arch are established. Since the blast load and soil 


load are expressed in terms of a general function ) Po sin = one set 


of moment equations suffices for both analogous loads. 
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3.5.1. General Equations for the Right Reactions of Circular-Segment Arches 
Since the circular-segment two-hinged arch, Figure 3.2, is one 


degree indeterminate, one equation based on Castigliano'’s theoren, 

Reference (4), is used for the determination of the reactions. The horizontal 
thrust at the right reaction is chosen as the redundant. It is assumed that 
the arch rib may’ be treated as a thin ring, making it unnecessary to consider 
the theory of curved beams. The effects of rib shortening also are not 
considered. 


The following geometrical relations are derived in accordance with 


Figure 5.2: 
Figure 3.2a BC = 2R sins? 
AB = BC cosa( 6 - ~) = QR sinsP cosa{ 6 - 9) 
AC = BC sins(@ - 9) = QR sinsP sins(@ - 9) 
BD = 2R sinsé@ 
Figure 3.2b DF = 2R sins(@ - 8B) 
EF = 2R sina{@ - B) cosa{6 - B) 
= moment arm of P about D 
Figure 3.2c FG = R sin (9 - 8) 
= moment arm of P about G 
Figure 3.ed AB = moment arm of V about C 


AC moment arm of H about C 


DB = moment arm of V about D 


In Figure 3.2d, the total strain energy U in the arch rib due to 
the moment M caused by the radial load P acting at a point on the rib 


located by B may be expressed as 





Li 


oe, ee 


Since the abutments are assumed not to move, the contribution of the right 


reaction H to the strain energy is zero and may be expressed as 


sre | BH Rae = 0 (a) 


Two equations are required to express the moment in the arch rib; 


% 


one equation for angle P less than or equal to f, and one for 9 equal to or 


greater than. For 0 <9 SB, 


M, = V(AB) - H(AC) 
= 2VR sin 39 cos3(@ - ) - 2HR sins? sins{ 6 - 9) () 
b 
3 
Sp. -OR sins? sins(@ - 9) (c) 
For B < ? < 6, 
M, = V(AB) - H(AC) - P(FG) (d) 


= 2VR sins? coss(@ - 9) - 2HR sind? sin3( 6-9) - PR sin(?#) 


aS 


sq = -eR sins? sins({ 0 - 9) (e) 


Substituting Equations (b) through (e) into Equation (a), 


o= fu, se naps |'m, ge Rae 


B 

- | [ 2VR sins? cos3{ 6-9) - 2HR sins? sinz{ 6-9) ] 

o 
[ -2R sinsp sins(@-9)] Rap + 


) 
[ [ 2vVR sins¢ cosa{ 6-—) - 2HR sins? sins 6-9) 
B 
- PR sin(9-8) J(-2R sinsp sins a) ] Rag (f) 


Solving Equation (f) for H yields 
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| “fo VR’ sin” 5 8 sin( SY =) cos( LY ap 


-2 i PR? sin(?-68) sin : sin(F)ap 


ee (g) 
uf a? ain? 2 esa & eed 


Referring to Figure 3.c, the sum of the moments equated about 


the left reaction gives 


2pR sin( SP) cos( 28) 


Ve P( EF) _ 2 ns sin( 6-8) (10) 
DB OR si @ - @ 
sin; 2 sins 


Substituting Equation (10) into Equation (g), 


are je-% yZ sine - s cos | + cos 5 t cos(# = 6) 


- F con? - p) + (25%) sin(Z - p) | | 


~ + 6 cos“(¥) - 4 sin - cos - 
11) 


Equations (10) and (11) define the right reaction of any circular-segment 


two-hinged arch acted upon by a radial load P located at any angle 6B. 


3.5.2. Moment Equations for 180 Degree Arches 


For a semicircular arch, the central angle 9g is ", and 


Equation (11) reduces to 


a °) 


— (sin B + (B - 2) cos 6] (12) 


Tt 


Using Equation (12), the horizontal right reaction may now be 


found for the harmonic load component p sin 
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HD 


For P = P} sin 4 Rd? = P) sin? Rd?, since @ = 4A, 
1 x pk 
H ai Pp) sin P(sin? + (9-x) cos~] Rag as 
For P = “Po sin am Rao 5 
1 x : : 2 
H = x | -p, sin2p[sinp? + (9-x) cos?) Rd? = = p.R 
x |. 2 al 
For P = oe sin 49 Rd@, 
1 K 
H=-— -p, sin39@[sin® + (9-x) cos?~] Rdg = 4 p-R 
Je 2 3 
The loading configurations Py, sin = are shown in Figure 35.1. 


It is also necessary to investigate the effect of the uniform 
antisymmetrical load P, and the uniform symmetrical load P, shown in 


Figure 5.1. For the uniform antisymmetrical load, 


x 
ae x 
H = < [? -p [sin? + (9-x) cos?] Ra? +e fs plsin? + (9-x")cosPJRaP = pR 
. 2 
For the uniform symmetrical load, 
x 2x 
_i| 35 dja 
H==5 -p[sin? + (9-2) cos] Rap + lhe pl{sin? + (9-x)cosP) Rap 
O — 
5) 


x 
a ae 3 
+ = te p [sin + (9-x) cosp] RaQ = 5 > PR 


2 


Equation (10) is used to find the right vertical reaction. 


With @ = x, Equation (10) reduces to V = 5p sin?. 





For P = 


For P = 


For P= 


%, 


For the 


Vs 


For the 
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il 


n 
P) sin = Rd?, V = af P) sin“? Rd? 
fo) 


ED IR 


x 
“Pro sin = Rd?, V = - | p, sine? sin? Rd? = O 
O 


6 


Ps § 
“P, sin _ Rd?, V = a P, sind? sin? Rd? = O 
oO 


uniform antisymmetrical load, 


* x 
2 | -p sing Rap + 3 | p sin? Rd? = O 
Oo pi 


uniform symmetrical load 
Ou 


A 
ved]? -p oin Rd +3 ],2 p sing Rap 
Oo = 


3 


514 
Qn p sin? Rd? = O 
z, 


i 
oe 


The value of zero for the last four loads was to be expected from 


the symmetry of the arch and the loading configurations. It is interesting 


to note that the two-hinged semicircular arch is statically determinate for 


the last four loadings. 


All the 


above values for the right reactions for a 180° arch are 


tabulated in Table II. 


The moment at any point onthe arch in Figure 3.1 located by ? 


due to a distributed load p, sin = RdB acting at B is 


dM = -p, sinB[R sin (9-6) ] Rap 


The total moment at 9 is 





p 
M = -H(AC) + V(AB) - [ p, sin8 (R sin(9-6)] Rdp 
O 


8 TR 
w-=E | on sin 8 oin( $8] +E | on sin $ cool 2 %) | 
* | & sine sin” B - conn (8. in ) |" 
O 


=) é p,R* K - (A-29) cos? - 3 sin | 


A similar procedure is used to find the moment at P due to the 
Other loading configurations in Figure 3.1. The results are tabulated in 


Table III. 


5-3-5- Moment Equations for 120 Degree Arches 
When the central angle of the arch is 120 degrees, Equation (11) 


reduces to 


{ED Z sin 3 - 3 cos 3 







B-=—> 
= - B) - bcos (x-6) +—— 


K 
+ COB = - : cos(= - 





4 ox a. x 
A i ae = - 4 sin = cos = 


3 3 3 


P | (comp +t oie 2- Bn + 3 -(2 8 oR cosp 


se sinp +p cos B - $B sin | 
NE 


re | 1.83979 - 2.70578 cosB + 0.16665 sing 


+ 1.59330 B cos B - 0.91990 B sine | (13) 





ce 
Equation (13) may be used to find the right horizontal reaction 
of an arch with a 120 degree central angle acted upon by a concentrated 


radial load P located at any angle 6B. The reaction for Py sin = Rd? is 


found as follows: 


Ox 
H = 3 P) sin - Dp | 1.83979 - 2.70578 cosP + 0.16665 sin? 
O 


x 


+ 1.59330 PcosP - 0.91990 pein | Rd? = 0.53076p,R 


The horizontal reactions for the remaining harmonic loads and the 
uniform loads are found in the same manner, and are set forth in Table II. 
Equation (10) is used to find the vertical reaction. 


For Py sin = Rd? , 


ox 
= jee ee 2 ap 
V = , z |», sin Zo sin (3 ~) Rap 0.6 p|R 


2 sin = 
D 
The remaining vertical reactions have been found, and are stated 
in Table II. 


For a segment RdB at B loaded with Py sin & Rd 68 , the moment 


at ~ in Figure 5.1 is 


dM = -p, sin 2 6(R sin(9-6)) Rap 
and the total moment at 9 is 


M 


-0.53076 p, Rfac] + 0.6p, R(AB] 


: 3 
- |p, sin 2 ptrsin(9-s)} Rap 
Oo 
=~ - 4 
= -0.53076 p,R E sin = sin( 2 —) | 


2x 
+ 0.6 p,R E sin - cos(2——) | 


‘® 
- | Pp, sin 2 B (R sin(9-6)) Rap 
= (0.78499 - 1.35965 sing - 0.78499 cosy + 0.8 sin _ p,R* 
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The moment equations for the remaining load configurations shown 
in Figure 3.1 for an arch with a 120 degree central angle have been found 


and are contained in Table III. 


3.3.4. Procedure in Determining Analogous Loads 
The moment equations tabulated in Table III have been evaluated 
for every ten fee The results are tabulated in Table IV and shown 
in Figures 4.5 and 3.4. These values are for static loads. They do not 
indicate the dynamic moment which may exist at any time in the arch. 
Inspection of Figures 3.3 and 3.4 shows that the effect of the 
analogous antisymmetrical load is almost identical to that of the load 
Po sin = be Po has a value of one, then the value of the antisymmetrical 
load must be slightly less than one if the effect of each is to be the same. 


However, the antisymumetrical load must also include the effect of the loads 


Py sin - and Ps sin a - Therefore, the moment that) P, sin a causes 


at a point 6/4 is first determined, and subsequently a value of the anti- 
symmetrical load is found which will cause the same moment at 6/4. The 
point 6/4 is selected since the unsymmetrical blast load produces the maximm 
effect at this point. The time must remain constant for each numerical 
evaluation of p(t), as described below. 

The actual blast load on the arch is p(t). At any given time 
the load profile p (9, t) can be determined and analyzed to give a distributed 


e] 
its effects may be found at 6/4. An analogous load can be found which 


load) p(t) sin > . By assuming that ) P( t) sin nt acts statically, 


causes the same effect. By repeating this operation at various times, the 
pressure-time diagram of p(t) can be determined. This load must then be 


applied dynamically to determine its actual effect on the arch. 





2k 
In determining the static load of the earth, the analogous load 
is found by equating the moments at 6/4 of the distributed static load 


6 
pressure is not a dynamic load, only one evaluation is needed. 


» PB, sin a with the moment of the antisymmetrical load. Since the earth 


Further inspection of Figures 3.3 and 3.4 shows that the effect of 
the distributed loaas ) De sin — on the central third of the arch may be 
represented by the symmetrical load which acts in compression on the central 
third of the arch and in tension on the outer thirds. This load is 
designated Pea If the distributed load is a time-dependent load, the same 
procedure is used in finding the pressure-time diagram for pC t) as that 


used for P( tle 





2D 
TABLE I 
COEFFICIENTS FOR RESOLVING VERTICAL AND HORIZONTAL PRESSURE INTO 


RADIAL PRESSURE 


Angle Radial pressure p(t) in terms of vertical pressure 


‘ Assumed ration p/P, 

One> 0.50 0.75 1.00 

O 0.250 0.500 0.750 1.00 
10 0.273 0.515 0.758 1.00 
20 0.338 0.559 0.779 1.00 
30 0.438 0.625 0.813 1.00 
LO 0.560 0.707 0.853 1.00 
50 0.690 0.794 0.897 1.00 
60 0.813 0.875 0.938 1.00 
70 0.913 0.942 0.971 1.00 
80 0.978 0.985 0.993 1.00 
90 1.000 1.000 1.000 1.00 
100 0.978 0.985 0.993 1.00 
110 0.913 0.942 0.971 1.00 
120 0.813 0.875 0.938 1.00 
130 0.690 0.794 0.897 1.00 
140 0.560 0.707 0.853 1.00 
150 0.438 0.625 0.813 108 
160 0.338 0.559 0.779 1.00 
170 0.273 0.515 0.758 1.00 
180 0.250 0.500 0.750 1.00 
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(bd) Uniform Symmetrical and Antisymmetrical Loadings 


Figure 3.1 Various Load Configurations 
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Figure 3.2 Geometry of a Circular Arch Segment 
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Figure 5.5 Moment vs Angle ? for 180 Degree Central Angle 
Two-Hinged Arch 
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Uniform Antisymmetrical 


pR° Moment, inch pounds 





Angle @, degrees 


Figure 3.4 Moment vs Angle 9? for 120 Degree Central Angle Two-Hinged Arch 
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CHAPTER IV 


BLAST PRESSURE LOADING 


4.1 Pressure-Time Relations at a Point 

As stated in the general discussion of the variables, Section 2.l, 
the blast pressure acting on the surface of the ground above the point of 
interest causes a vertical pressure at that point. Accompanying this 
vertical pressure is a horizontal pressure of varying magnitude, depending 
on the soil. The horizontal pressure may be as low as fifteen percent in 
loose, dry soil, fifty percent in moist soil, and one hundred percent in 
saturated soil, Reference (1). One of the purposes of this thesis is to 
study the effects of this variation. 

The vertical pressure is propagated downward with a velocity 
equal to the seismic velocity of the soil. At the point of interest, the 
time for the pressure to increase to a maximum is a function of the vertical 
transit line. The effect of varying the rise time will be studied. 

It has been shown in Section 3.1 that the horizontal pressure 
may be expressed in terms of the vertical pressure. Vertical pressure-time 
relations at a point on the arch in terms of surface pressure are there- 
fore required to define the total radial load profile PA, t) acting on 
the arch. The equations below are derived from the geometry of Figure 4.1. 

The velocity of the wave propagation is: 


C_U 
8 


v= ore (14) 


where Vv velocity of propagation of the earth pressure wave 


C = seismic velocity of the soil 


U = surface blast wave velocity 





3h 
The time required for the pressure wave to arrive at any point 


onthe arch rib is: 


t _ R - R sin(@ - 8) (15) 


a Vv 
where @? = angle locating a point on the arch rib 


R = radius of the arch 


o = time measured after initial contact of the pressure 
wave with the arch rib 
C., 
B = arc tan a 


The time required for the pressure at P to rise to its peak 


is; 
R - R sin? + i. 
t= ne » fort = t, (16) 
R - R sin? + #H ‘ 
eC acm fort.=2t, 
where H. = depth of cover on the crown 
c. = rise time of the pressure 
, t, = vertical transit time of the pressure wave to ? 
R - R sin? +H 
t, = : 
t C 


s 
The vertical pressure-time variation at a point on the arch rib 
is shown qualitatively in Figure 2.2. The pressure-time diagram is assumed 
to be a straight line in the rise part of the diagram and to have an ex- 
ponential decay rate after reaching its peak value Pa The total impulse 
of the force pulse at any point must be the same as the total impulse of 
the surface force pulse. This means that as the rise time increases, the 


decay rate must also increase in order to maintain the same impulse in 





bie. 


the force pulse. It is assumed in this study that the decay rate at the 
shallow points in the soil considered in this thesis is the same as the 
decay rate of the surface force pulse. This is a reasonable assumption 
since the maximum unsymmetrical loading occurs within the first few milli- 
seconds after the blast wave begins to act on the arch. 

The decay rate for a side-on overpressure Pp, baving a positive 


phase duration t, may be obtained from Reference (6). The decay rates for 


d 
the two overpressure levels studied in this thesis are shown in Figure 4.2. 


For the 100 psi overpressure, t., is assumed to be 1500 milliseconds, while 


d 
for 200 psi overpressure, t, is assumed to be 1640 milliseconds. The 


d 
weapon yield is 1 MI. 

The complete vertical pressure-time diagram for any point can be 
determined only by numerical evaluation of the pressure on that point at 
various times. A sample calculation of a pressure-time diagram follows. 

Assume that a 180 degree central angle arch with a radius of 
twelve and one-half feet is subjected to a 200 psi overpressure blast wave 
originating from a 1 MI weapon. The duration of the positive phase of the 
loading is 1640 milliseconds. The arch is buried in soil having a seismic 
velocity of 1000 fps, with zero depth of cover on the crown. The time of 
rise of the earth pressure at a point in the soil is assumed to equal the 
vertical transit time of the earth pressure wave. The ratio of the 
horizontal to vertical pressure is assumed to be 0.25. 

Reference (7) gives the following equation for the shock velocity 


U of the surface blast wave: 


Use (1 + 6p/TP,)* 


where Cae ambient sound velocity ahead of the shock 





36 


peak overpressure behind the shock 


~~ 
iH i 


ambient pressure ahead of the shock 


ee 
~°. U = 1130 (1 + 1200/103)* = 4023 fps 


es 


me 


UC 
By Equation (14), v= ——=—=—===— 

\°. Ap \| 10007 x 4023" 
By Equation (15), 6B = arctan (1000/4023) = 13°58! 


t= ora (2 - sin(? - 13°58')] , seconds 


= 12.8 (1 - sin(9 - 13°58'] , milliseconds 


= 2d : 
By Equation (16) t. 7300 | 1 - sin?) , seconds 
= 2.5 (1 - sin®?), milliseconds 
Let 9 = 40°. ‘Then 


t = 12.86 (1 - sin(40° - 13°58')] = 7.2 ms 


& 


t= 12.5 (1 - sin 40°) = 4.5 ms 


Therefore at = 40 degrees, the earth pressure wave arrives 
7.2 ms after the wave initially contacts the arch rib at 9 = 90° + 13°58'. 
The vertical pressure rises at a rate of 200/4.5 = 45 psi per ms. At 8 ms, 
the pressure is (8 - 7.2)45 = 36 psi. At 10 ms, the pressure is 
(10 - 7.2)45 = 126 psi. At 11.7 ms, the pressure reaches a maximum value 
equal to the surface overpressure, or 200 psi. Thereafter the pressure 
decays at a rate assumed to be equal to the decay rate for the surface 
overpressure as shown in Figure 4.2. Inspection of Figure 4.2 shows that 
the decay may be represented accurately by a straight line down to a value 


of 4 P/OP = 0.65 and dt/t, = ~ 0.02. The pressure decrement is then 
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0.35(200) = 70 psi, and the time interval is 0.02(1640) = ~ 35ms. So the 
pressure decays at a rate of 70/35 = 2 psi per ms. Therefore at 12 ms, the 
pressure is 200 - 2(12 - 11.7) = 199 psi. At 14 ms, the pressure is 

(199 - 4) = 195 psi. This procedure is continued until a time greater than 
(35 + 11) ms is reached. Then the time interval is expressed as a fraction 
of t, and AP/AP “is read from Figure 4.2. The complete pressure-time 
relation for ? = 40 degrees is recorded in Appendix B, Table B5. 

To obtain the complete vertical pressure-time relations for the 
entire arch, the above procedure is accomplished for every ten degrees on 
the arch. Table B5, Appendix B, shows the vertical pressure-time variation 
for the entire arch up to 200 milliseconds. The unsymmetrical loading is 
readily seen. For example, at ten milliseconds the windward side of the 
arch has pressure acting on its entirety, whereas the leeward side is not 
fully loaded until after sixteen milliseconds. The unsymmetrical loading 
shifts from the windward side to the leeward side during the time the blast 
pressure acts. This effect may be seen at thirty milliseconds. The 
pressure on the windward side, greater than ninety degrees, averages some 
eight pounds less than the pressure on the leeward side. The effect of 
the decay rate may be seen by comparing the pressure at the crown with 
that at the abutments. For example, at twenty-four milliseconds the pressure 
at the crown has reached its maximum and decayed to approximately seventy- 
five percent of the maximum surface pressure, while the pressure at the 


abutment is approximately equal to the maximum surface pressure. 


4.2. Harmonic Analysis of Radial Blast Load with Time a Constant 


Based on the vertical pressures tabulated in Table B5, a radial 


load profile p,.(?, t) can be determined for any instant of time given in 
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that table. After p_.( 9, t) has been defined, the distributed load 


8 
seconds, is given in Tables V and VI. 


» p(t) sin ae? can be evaluated. A sample calculation, for t = 10 milli- 


In Table V, the data in column "p(t) " is obtained from Table B5, 
while that in colum "p_ coefficient" is derived from Table I for p/p, = 
0.25. The data in colum "p_(P, t)" of Table V is derived by miltiplying 
p( t) by the p_ coefficient. The resulting radial pressure load profile 
pl 9, t) is shown graphically in Figure 4.3. The coefficients p,,( t) of 
Equation (9) are determined by numerical integration. This is done by 

nm 


multiplying the load ordinate at each ten degree interval by S sin Ta? 


adding the results and miltiplying the sum by > - The resulting dis- 
tributed radial loading ) p(t) sin =? , with t = 10 milliseconds, is 
shown as a dotted line in Figure 4.3. 

Table VI contains the harmonic analysis of the same blast wave 
acting on a 150 inch radius arch with a central angle of 120 degrees. The 
radial pressure p.( 9, t) in Table V from 30 degrees through 150 degrees is 
used in Table VI from O degrees through 120 degrees. The resulting dis- 
tributed redial loading p(t) sin a is shown in Figure 4.3. 

Inspection of the data in colum "p_.( 9,t)", Table VI, reveals that 
there is a radial pressure of ten or more pounds acting at all points on 
the arch rib. Therefore, a uniform all-around pressure of ten pounds may 
be subtracted from the radial pressure before the harmonic analysis is made, 
since it is desired to resolve harmonically only the radial pressure causing 
flexure in the arch rib. The colum ape adjusted” contains the remaining 
redial pressure. 

The coefficients p(t) of Equation (9) are found by numerical 


integration. This is done by miltiplying the load ordinate at each ten 


. TR 
degree interval by is sin eae » adding the results, and multiplying by = ; 





dY 


The resulting distributed load, p(t) ) p(t) sin oe » with t = 10 milli- 


seconds, is shown as a dotted line in Figure 4.3. 
It is assumed that the first three terms of the series 


e) 
of this thesis. 


» Pp, sin uz adequately represent the loading on the arch for the purposes 


The radial load on the 180 degree arch at a time t = 10 milli- 
seconds is; 


p(t) = 166 sin? - 15 sin29 - 26 sin3@ 


Since the time is constant, the load may be considered to be a static load. 
An analogous load may be determined which will have the same effect on a 
given point of the arch. Let the point be located at P = 45 degrees. 


e) 
negative moment at this point. Using the coefficients in Table IV, the 


inspection of Figure 4.3 shows that the components) Po sin ax all cause 


moment of the load is 


M= [166(-0.02262) + 13(-0.33333) + 26(-0.08837) ] Ro 
ai[23. 75s Nekoi” = -10.562— 
A uniformly distributed antisymmetrical load acting in on the 
windward half of the arch and out on the leeward half produces a moment of 
o.4ikeep RO at P= 45 degrees. Therefore 


2 
10.388 = 25 psi for a time t = 10 ms 
O-41420R 


a 
However, the actual load on the arch is time-dependent and consequently De 
must be a time-dependent load, p( t). To define the pressure-time diagram 
for p(t) » the load P.(, t) on the arch is harmonically analyzed at various 
times. 


The analogous load for the crown may be found in a similar manner. 


The moment caused by the above load at the crown is 





hO 


M s_ [166(0.03540) + 26(0.125)] Ro 


= 9.15 ae inch pounds 


A symmetrical load P, acting in compression on the central third 
of the arch and in tension on the outer thirds of the arch produces a 


moment of 0.1547 pR at the crown 


% 


2 
ne p(t) c= ae = 59 psi, t = 10 ms 
0.1547R 


4.3 Blast Lo Approximations 

Appendix B contains the results of numerous harmonic analyses of 
the blast load on an arch arising from a 1 MI nuclear weapon. The analyses 
were made and analogous loads found as described in the preceding section. 
These analyses were made in order to determine the effects of the variables 
discussed in Chapter II on the loading the arch receives, and to provide a 
basis for idealizing the blast load on the structure without the necessity 
for actually harmonically analyzing the blast load. 

The blast load mist be defined in relation to the effect being 
studied. In this thesis, the effect studied is the deflection of the arch 
due to flexural stresses, the area of primary interest being the haunch of 
the arch. The effect of the unsymmetrical blast load on the haunch of the 
arch if most conveniently related to an analogous uniformly-distributed 
load which acts inwardly on one half of the arch and outwardly on the other 
half. The effect of the blast load on the crown is conveniently related 
to a uniform load acting in compression on the central third of the arch 
and in tension on the outer thirds. Values of such loads for various 


assumed conditions are shown in Tables Bll and Bil2. 
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After the blast load has been related to the specific effect 
being studied, a load-time curve must be established to show the variation 
of the load with time. It will be noted in Tables Bll and Bl2 that the 
analogous loads are for a specified time only, and that they vary appreciably 
with time. Appendix B contains a few typical analogous load-time curves. 
Similar plots were made of all analogous loads. 

Study of these load-time curves revealed that the peak values of 
the analogous loads are some fairly uniform function of the maximum surface 
overpressure Pi for all assumed values of the variables. Im addition, the 
shapes of these load-time curves are approximately initially peaked tri- 
angular load pulses, which decay rapidly to a relatively small value and 
then continue to decay at a much slower rate. The duration of the rapid 
decay phase of the loading corresponds approximately to the time required 
for the pressure to engulf the arch. The duration of the slow decay phase 
of the loading is relatively much longer. 

The peak values of the analogous loads found in Appendix B are 
shown in Table VII and plotted in Figures 4.4 and 4.5. The data in 
Appendix B and Figures 4.4 and 4.5 may be qualitatively summarized as 
follows: 

(a) For any given overpressure, the peak value of the 

analogous load for the haunch increases as the value of 

the seismic velocity of the soil increases. This is ex- 

plained by the fact that a higher seismic velocity allows 

the pressure to penetrate the soil to a greater extent 

and thus a larger segment of the arch on the windward side 


is loaded before the leeward side is loaded. 





(bd) 


(a) 


oe 


For any given overpressure, the peak value of the 
analogous load for the crown decreases as the seismic 
velocity of the soil increases. This is explained by 
considering the effect of the load profile P_.( 9 t) on 
the integral ie p_.(?, t) sin = Rap - When C. has 
a relatively high value in comparison with U, a greater 
portion of the segment from ?) to Po is located on the 
windward side of the arch and the values of the sine 
function of the angle ? are both positive and negative. 
When C. has a lower value in comparison with U, the 
load profile p,.(?, t) is located more nearly in the central 
third of the arch, where the value of the sine function 
of the angle 7? is negative. This increases the numeri- 
cal value of the analogous crown load. 

The higher ratios of PB, /Py increase the peak value 
of the analogous load for the haunch. This is explained 
by the fact that the area of the total radial load profile 
p (9, t) is greater for the higher ratios, and thus the 
total unsymmetrical load is greater. 

The p/P, ratio has little effect on the peak value 
of the analogous crown load. The peak value of the crown 
load is primarily a function of the location of the load, 
i.e., the central third of the arch. In this area, the 
value of the coefficient for determining the resultant 
of the horizontal and vertical forces is approximately 
one for all values of P,/ Pays Consequently, changing the 


ratio has little effect on the radial load. 
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(e) The assumed rise time has relatively little effect 
on the peak value od the analogous loads. The peak 
values of both loads are reached very early in the load- 
ing of the arch when the pressure wave has not penetrated 
the soil to any extent. For shallow depths, the rise 
time is" practically equal tothe rise time of the surface 
pressure, which is assumed to be instantaneous. 

(f) Although the higher ratios of P,/P, cause a higher 
peak value of the analogous haunch load, they actually 
cause a less severe dynamic load effect. If the ratio 
is high, the pressure, after the pressure wave has en- 
gulfed the arch, is more nearly a uniform compression, 
causing less flexural action in the arch. On the other 
hand, if the ratio is low, the uniformity of the compres- 
sion is reduced, increasing the flexural action in the 
arch. For the higher ratios of p/P, it is found that 
the analogous loads change direction, thereby creating 
a negative phase of loading. The peak value of the 
pressure in this negative phase is small, but could pos- 
sibly cause an outward movement of the crown of the arch 


since it has a relatively long duration. 


4.3.1 180 Degree Arches 
Since the blast loading may act over all, or over part, of the 


arch segment, one approximation to the peak analogous loads may be obtained 
by considering the blast loading acting on the arch in such a way that the 


maximum value is given to each component of loading in the series 





yy 


3 pect) sin a - The ratio of the horizontal pressure to the vertical 
n 


pressure which causes the most severe dynamic loading is 0.25, therefore 
this value is used in the following derivations. The effect of the rise 
time of the pressure and the time of arrival of the pressure wave is not 
considered. The vertical pressure at all points on the arch is assumed to 


be equal to P . The unit radial pressure at any point is, by Equation (3), 
p(t) = (0.25 + 0.75 sin“) P_ 


Using Equation (9), 


— ne 
p(t) oR | P,.(? t) sin 5 Rd 


maximum P, ( t) 


2P a D) 
=| (0.25 + 0.75 sin*9) sing RAP = 0.955 PL 
.@) 


oP PA 
maximun Pp(t) = ae (0.25 + 0.75 sin“) sine? RdP = -0.398 Pa 
2 


On 
op = 
maxi nun pz(t) = TE NS (0.25 + 0.75 sin“) sin3? rdP? = -0.408 Sa 


5) 

The maximum effect of the blast loading which may occur, under an assumed 
ratio of P,,/Py = 0.25, can never exceed the effect of the sum of maximum 
values of the harmonic components found in this manner. The value of the 
individual maximum values found, and the effect of the sum, are always 
reduced when the seismic velocity of the soil, the shock velocity of the 
surface blast wave, and the pressure-time relationships are cansidered. 

If these components acted statically and separately, they would 


produce the following moments at the leeward haunch and crown: 
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M = - | 0.025(0.955P,)R* +0.333(0.398P_)R° + 0 .088(0.4OSP, JR? 


2 
= -0.190 PR 


2 2 2 
B= 0.035(0.955P_)R + 0.125(0.408P_)R = 0.085 PR 


The peak value of the analogous haunch load is 


2 
0.190 POR 


p(t) = ———— = 0O.46P 
2 O.41h Ro = 


The peak value of the analogous crown load is 

0.084 P RS 
p(t) = a = O.5% P 
These values of the analogous loads can never occur, and, for 

want of a better name, they are designated as the maximax analogous loads. 

They are shown in Figure 4.4 so that they may be compared with the analogous 

loads found by actual harmonic analysis. Such comparisons indicates that 

these approximations are conservative. The following reductions are 


arbitrarily made: 


Maximum p(t) = 0.37 P. (17) 
Maximum p(t) = 0.45 ise (18) 


Equations (17) and (18) are plotted in Figure 4.4. 

To obtain an idea of the magnitude of these analogous loads, they 
may be related to a uniformly distributed load acting on a simply supported 
beam. Inspection of Figure 3.3 shows that the moment caused by the analogous 
load P, on the haunch of the arch is very similar to the simple beam 


moment. The moment at 45 degrees caused by the analogous load is 


M = 0.37(0.41422) PR 
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Assume that the simple beam has a length equal to the length of half the 
arch, or L = - . The moment is wL"/8. Equating this to the moment in 


the arch, it is found that: 


8 x 9.37 x 0.41433 PRS 


2 
(=) 


= 0.5 Pa 


A similar analogy may be drawn for the analogous load acting on 
the crown, with the exception that the length of the simple beam is now 
one-third the length of the arch, or = - The moment in the arch at 


b) 
90 degrees due to the analogous load is 


cq 
tt 


0.45(0.1547) PRS 


8 x 0.45 x (0.1547) PR 


z 
x 
“wolf 


(By 
p) 
= 0.5 Ps 
The above maximax loads occur within approximately two to four 
milliseconds after the loading begins to act on the arch. This fact is 
evident by inspection of the data in Appendix B. These loads decrease 
rapidly since Pp(t) becomes practically zero when the pressure wave has 
engulfed the arch, and pz(t) also decreases. However, p,(t) and p(t) are 
still acting. It is now assumed that the pressure has completely engulfed 
the arch. Since the pressure has engulfed the arch, the all-around uniform 


pressure may be subtracted from the expression for p_.(9, t) « Then 


p(t) = 0.75 sin'9 Pp. 
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This radial pressure may be resolved into harmonic components. It is 


found that 


2P x 
ml 2 ; ” 
p, () me. J (0.75 sin ?) sin? Rd? = 0.637 Ee 


oP . 2 
p,(t) = xR i (0.75 sin ?) sin3? Ra? = -0.127 e 


The moment in the haunch at 45 degrees caused by these two 


components is 


7 == | 0.003(0.657P,)R* + 0 .088(0.127P, )R° | = -0.026 PR 


The haunch analogous load is 


p(t) = sayy «= (0.06 PL (19) 


The components p(t) and p(t) act for the total duration of the 
surface blast pressure. Reference (2) states that the duration of the 


positive phase of the side-on overpressure may be considered as 


2/3 
t, = 0.40 (200) 5, 1/3 (20) 
Pp 
m 
where W = weapon size in MT 
t, = the effective duration in seconds of a triangular pulse 


having the same impulse as the actual shock. 
The total duration of p, (+) and p(t) is considered to be t,. 


The time required for p,(t) to become practically zero is found 


by using Equations (15) and (16), with the angle 9 = 0. 
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This time is defined as: 


+ q 2 te + t.. = duration of initial triangular impulse 


oo R +) Rein pb pga - 2 
= CU oC" » with be 2 t, 


The time t! defines t for the load p,(t) = 0.06p_ 


4.3.2 _120 Degree Arches 
Approximations of the blast load acting on a 120 degree arch may 


be derived in the same manner as that used in finding the approximations 
for the 180 degree arch, as described in Section 4.3.1. It is necessary to 
modify Equation (3) by substituting (9 + 30) for 9, since in Equation (3) 
? is measured from a horizontal radial line joining the arch abutment with 


the center of the circular segment. Equation (4) becomes 


[0.25 + 0.75 sin“(9 + 30)] Py 


It + - sin“ + 33 sin? cos? + = cos"9 Pia 


Let p(t) = | PA? t) sin oe Rd 


p.(9, t) 


Then the maximum possible values of p(t) are 


Cn 


j2 ( t + = sin°9 + BP) sin? cos? 


“Ble 


Maximum p, (t) 


+ = cos9) sin 2 P Ra? 


= 1.103 Pa 





2) 


On 
op 3 
Maximum Pp(t) = = r d ( ; 7 = sin“ + 3.5 sin? cos? 
= K 
~ B 
+ & cos“¢ ) 8in49 Ra? 
= -0.505 sae 
*, hn 
oP fe NE 
_ on [92 oy ae 
Maxi mum p(t) a i ( me = sin P + 3 sin? cos? 
? 9 
+ - cos“ ) sin 4 Q RaQ 
= -0.417 P 


The harmonic components p(t) sin = would produce moment at the 
haunch and crown as shown below, if the blast load is considered to act 
statically. The moment coefficients are obtained from Table IV, for P = 30 


and 60 degrees. 


M = - | 0.009(1.103P, )x® + 0.125(0.505P_)R© + 0.050(0 417)? 
= - 0.094 PRS 

M, = 0.014(1.103P_)R° + 0.051(0.417P_)R* 
= 0.037 PR 


Using the moment coefficient for ~ = 30 and 60 degrees in 


Table IV, the analogous loads are found to be 


0.094 PR 


Maximum p (t) 
= 0.1547R° 


0.61 P 
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2 
0.037 Pa 


0.06556R° 


Maxi mum p(t) 


0.56 Ea 


These values are plotted in Figure 4.5 so that they may be 
compared with the “actual analogous loads found by harmonic analysis, as 
given in Appendix B. Since they are conservative, the following arbitrary 
reductions are made; 


Macc mum p(t) = 0.45 Ee (21) 


Maximum p(t) = 0.41 Pa (22) 


These equations are plotted in Figure 4.5. 

These loads may be related to a simple beam load of 0.5 Pa by the 
same method as that described in section 4.3.1 for a 180 degree arch. 

After the Suey rive wave has engulfed the arch, there is a uniform 
all-around pressure acting which must be subtracted from the radial 
pressure in order to determine the loading which causes flexural action. 


This uniform pressure is 


— 2 cs 
p= (0.25 + 0.75 sin 30) Pa = - Pa 


The remaining radial pressure which is to be resolved harmonically is 


Sr 
p,.( 9) - P= ( ' + = sin’ * 2x sin? cos? + -- cos“9)P_, < - ee 


( = ra ~ - sino + 33 sin? cosP + 2 cos“) Pa 


Then 





ppl 


On 
oP =< 
p, (+) = a J y (- - + = sin“? + 5% Bin? cosP + = cos“9) sing? Rd 
5, 
= 0.546 P 


On 
OP = 
p(t) = a i (- - + _ sin-@ + 3.5 sin? cos? + = cos?) sind Rd? 


= -0.020 P 
ma 


These loads would produce the following moment at the haunch and 


crown of the arch if they acted statically: 


2 2 
w=: | 0.009(0.5K6P,)R + 0.050(0.020P_)R 
= -0.006 P_R* 
m 
Mr = 0.014(0.546P_)R® + 0.051(0.020P_)R® 
2 
- 0.009 PER 


The analogous antisymmetrical load for the haunch is 

0.006 PR 

2) = a 

O.155R 

= 0.04 P 
m 


The analogous load for the crown is 


0.009P_R* 


ee) ee 
: 0.06556R 


= O.14 P 
m 


(23) 


(24) 
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The time t that defines the analogous loads in Equations (23) 
and (24) is the time required for the pressure wave to engulf the arch and 
rise to its maximum value at the abutments. This time is expressed by 
Equations (15) and (16). The angle 9 in both of these equations is thirty 


degrees. 


4.5.5 Summary 

The blast load acting on the arch may be defined by numerically 
evaluating at specific times the coefficient p(t) for the distributed load 
). p(t) sin #2 - The load-time curve for the analogous loads is obtained 
by equating the static effects of the distributed Load ) p(t) sin — to 
those of the analogous loads. This is the procedure employed in Appendix B. 
The procedure is extremely tedious, and, in fact, unnecessary, since the 
blast load is not accurately known. Therefore, a more convenient method 
such as that established in Subsections 4.3.1 and 4.35.2 is desirable. 

In these subsections, the blast load acting on the arch has been 
idealized to two initially peaked triangular force pulses. For investigat- 
ing the effect of the blast load on the haunch of a 180 degree arch, the 
peak value of the first pulse is 0.57? and the peak value of the second 
pulse is 0.06P _. For similar investigation of a 120 degree arch the peak 
value of the first pulse is 0.45 Py whereas the peak value of the second 
pulse is 0.04 P; in each instance the duration of the first pulse is the 
time required for the pressure wave to engulf the arch, and the duration 


of the second pulse is equal to the effective duration of the blast loading. 
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Figure 4.1 Geometry of Blast Loading on a 180 Degree Central Angle 
Circular Segment Arch 
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Figure 4.2 Overpressure as a Function of Time 
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Figure 4.4 Peak Values of Analogous Loads: Antisymmetrical p(t) and 
Symmetrical p(t) for a 180 Central Angle Arch 
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Figure 4.5 Peak Values of Analogous Loads: Antisymmetrical p(t) and 
Symmetrical p(t) for a 120 Degree Central Angle Arch 
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CHAPTER V 


PASSIVE EARTH PRESSURE LOADING 


If there were no soil surrounding the arch, it would deflect into 
the dotted configuration shown in Figure 5.1 under the unsymmetrical blast 
loading. However, the soil restricts this deflection by passive earth 
pressure and the inertia of the soil mass. To evaluate the passive earth 
pressure, the segment of the arch on the leeward side from the springing line 
to the crown is assumed to act as a retaining wall. The soil volume and 
the relative plane of slip are shown in Figure 5.1. 

Under actual conditions, the body of soil located on the windward 
Side of the arch would most certainly be affected by the blast pressure and 
would have to move downward and laterally into the arch to load it. The 
energy required to accelerate this mass would lower the energy level in the 
blast pressure. In addition, after the surface pressure wave passes the 
arch, there would be a vertical pressure of varying intensity acting downward 
on the soil body on the leeward side of the arch. This would increase the 
shear strength of the soil, the passive earth pressure, and the inertia 
force. Evaluation of these later effects would be extremely difficult and 


is omitted in this study. This results in a conservative approach. 


5.1 Unit Passive Earth Pressure 
- Rankine's formula for the total horizontal passive earth pressure 


force = acting on the vertical plane AV with a height h in Figure 5.1 is 


p = 4wn'k 


Pp 
1 + sin 9' ; 
where K = 1- eine” ~' being the angle of internal 
friction of the soil 
and w = unit weight of the soil. 
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The unit horizontal pressure along the plane AV in Figure 5.1 
varies directly as y varies, y being the distance below the ground surface. 
This pressure then may be expressed as 


P, = wyK | (25) 


51.1 180 Degree Arches 
For a semicircular arch with no cover on the crown and 9' = 30 


degrees, since y = R({1l - sin®), Equation (25) becomes 


p, = wR(1 - sing) 32 = 3wk(1 - sing) (26) 


The angle 9 is measured from the horizontal radial line which joins the 
abutment with the center of the circular segment. 

In addition to the horizontal passive pressure, it is necessary 
to consider the weight of the soil above the arch rib. This soil exerts a 


unit pressure 
Pp, = wR(1 - sin?) (27) 


Using Equations (1), (26), and (27), the resultant radial earth pressure is 


p(?) = 3(p, +P.) + 2 (Pp, - py) cos 20 
= 4 [wR (1 - sin?) + 3wR (1 - sing)] 
+ 5 [3wR (1 - sin?) - wR (1 - sin?)] cos 29 


QwR (1 - sin?) + wR (1 - sin?) cos 29 (28) 


In Figure 5.1, the angle B, which defines the inclination of the 


plane of shear AS with the vertical plane AV, is 


3 (90° + ~'), Q' being the angle of internal friction 
of the soil 


p 


4 (90° + 30°) = 60° for ~' = 30° 
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From the geometry of Figure 5.1, the total weight of the soil 


associated with the passive earth pressure is 
. 60) oe 8 een? (29) 
We= | 2 R(R tan 60) + -—- | w= 1. 


Equation (28) is for an arch with no cover on the crown. If 
there is a depth of cover A. on the crown, a uniformly distributed hori- 
zontal pressure Su, as derived from Equation (25), and a uniformly dis- 
tributed vertical pressure WH» must be added to the radial pressure acting 
on the arch. For an arch with a central angle of 180 degrees, the resulting 


radial pressure due to H,, as obtained by Equation (1), is 


i 1 : 
p,.(9) = 2 (Pp, + P,) + 2 (Pp, P.) cos 2D 


7 sh ? 
& (wi, + 3wi) + 3 (Swi, vii) cos 29 


eo wH + WH. cos 29 (30) 


The total radial pressure on a 180 degree arch with a depth of cover H. 


is found by adding the pressures derived from Equations (28) and (30). 


p.(9) = OwR (1 - sin?) + wR (1 - sin?) cos 29 


+ CwH + WH cos 2 
c c 


9-1.2 120 Degree Arches 
Equation (28) is valid only for a semi-circular arch with the 


angle P measured from the horizontal radial line. For an arch with no 
cover on the crown and a central angle of 120 degrees, Equation (28) is 
modified by substituting (P + 30°) for 9. This yields 


p.( 9) = QwR [1 - sin (9 + 30)) + wR [1 - sin(® + 30)] cos 2(9 + 30) 
= wR ( 2- 343 sin? - £ cos? - > cos 39 


2 - sin 2D +3 cos 29) (31) 
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If there is a depth of cover H. on the crown of the 120 degree arch, the 
resulting radial pressure due solely to i. is 
1 i 
p(%) = 2(p, +p) + 3 (p, - py) cos 29 


= 3 (3wH + wi.) + 2 (Swi, - wi) cos 2(@ + 30) 


Q 


Ske 


= 2 WH + 3 wi. cos 2® - WH. sin 2p (42) 


The angle ~ is measured from the radial line joining the abutment with the 
center. 

The total unit radial pressure of a 120 degree arch having a 
depth of cover H, on the crown is found by adding Equations (31) and (32). 
The result is 


PA?) = wR ( 2 - 223 sing - $ cose - $ cos > 
~ B ein 2p + $ cos 29) + 2 WH. 


+ 3 WH, cos 20 ~ oF WH sin 29 (33) 


The total weight of the soil associated with this force is 


W 
S 


2 
2 (R +H)? tan 69 - SE | Ww 


2 


0.542 wR” + 1.732 WRH, + 0.866 vale 


5.2 Distributed Earth Pressure Load 
In considering the unsymmetrical blast loading on the arch rib, 


it is advantageous to resolve this loading into harmonic components since 


the resulting distributed loads ) p(t) sin — are easier to treat 


analytically than the actual irregular blast load. This advantage does not 


exist for the passive earth pressure loading, since original analytical 
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expressions are available which may be studied in determining the effect 
of the passive earth pressure on the arch. However, resolution of the 
passive earth pressure into harmonic loadings does facilitate visualization 
of the fact that the earth pressure may be thought of as a distributed 
resistance ) Pp sin “= in the arch rib. In addition, a desirable wni- 
formity of approach is provided for the entire problem. 

Equation (28) expresses the unit radial pressure acting on a 
180 degree arch with zero depth of cover on the crown. The total force 
acting on the arch rib may be resolved into harmonic components by 


Equation (9). 


p_.(9) = QwR(1 - sin?) + wR(1 - sin?) cos 297 


See ~ 
Po | Se | p(%) sin == Rap 
x 
neu HD 
P, eS  Cewa(a - sin?) + wR(1 - sin?) cos 2P)sin a Ra 
= 
_2 [2 25? pap - 
Ss ape : P_.(?) sin Rd? = 0.509 wR 
x 
. 2. | 2 on? pao = 
Ps, aR J. P(?) sin <= Rap = 1.180 wR 


Equation (33) expresses the radial pressure acting on a 120 degree 
arch having a depth of cover H. on the crown. The total radial force re- 
sulting from this pressure may be resolved into harmonic loading components 


as follows: 





p,(9) = wn (2 - 252 ing - 2 cose - cos 4— 
- % ain 29 + 2 cos 29) 
+ wi (2 + = cos 2 - XB sin 29) 
a in BMD 
BP, “Fe J RCO sin —>~ Ra? 
P) ve [fF (2 - 422 353 cino - 2 cosp - 3 cos 39 
te 
V3 l 3 
= = sin 20 + 5 cos 29] sin SP Ra? 
vutt, [> (2 +3 con 20 - YB sin 29) sin 39 Rae | 
° 
= 0.111 wR + 0.864 wi 
La 
Pr = ss | vt) sin 39 RAM = 0.181 wR + 0.987 wi 
=FR ve) 2 s 
3 
Ae! 
Pp ne }? p.(9) sin Zp Rd? = 0.198 wR + 0.557 wi 
? 2m Yo * E = 
2) 


>-) Analogous Earth Pressure Load 
In order that the effect of the unsymmetrical blast load on the 


arch rib may be determined conveniently, the blast load is represented as 
ean analogous uniformly distributed radial load acting in compression on the 
Windward side of the arch and in tension on the leeward side. This greatly 


facilitates the study of the dynamic effects of this load since the problem 
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of dealing with a time-dependent unsymmetrical load is reduced to that of 
dealing with a time-dependent uniform load. The problem may be simplified 
further by letting this uniform load define the parameters of the dynamic 
load acting on the single-degree-of-freedom system. 

It is logical to extend the concept of analogous loads to the 


‘unsymmetrical passive earth pressure acting on the arch. In the first place, 
a desirable Berna ty of approach to the entire problem is obtained. In 
addition, by defining the earth pressure load as an antisymmetrical radial 
load acting in compression on the leeward side of the arch and in tension 
on the windward side, the value of this pressure can be considered as inte- 
gral part of the resistance of the arch rib, since the direction of the 
earth pressure is directly opposite to that of the blast pressure. 
Furthermore, this resistance provided to the arch by the soil may be used 

to define the resistance of the single-degree-of-freedom system. 

The loading components ) Po Sin t defined in Section 5.2 , 
coupled with the moment coefficients from Table IV, are used to define the 
analogous loads for 180 and 120 degree arches. The resulting analogous loads 
define the maximum, or yield, resistance given the arch rib by the soil, 
Since the passive earth pressure is evaluated at its peak value. The 
resistance is evaluated as a function of the moment in the arch rib ata 
point located by 6/4, the same point used to define the analogous load for 
the blast pressure. 

For a 180 degree arch with no cover on the crown, the resisting 


moment due to the soil is 
M= -0.023(0.948 wR)RO + 0.333(0.509 wR)R° + 0.088( 1.180 wR)R® 


= 0.252 wR? 





The analogous load is 


«Oe wR? 


= O.4lh Ro 


0.611 wR (34) 
For a 120 degree arch with a depth of cover Be on the crown, the 


resisting moment due to the soil pressure is 


M = -0.009(0.111wR + 0.864wH )R° + 0.125(0.181WR + 0.987wH, )Ro 


+ 0.050(0.198wR + 0.557WH, )R® 


2 


0.034uR? + O.1uliwH R 


The analogous load is 


O.034wR? + O.144wH 


p ama STE ST aT ar 
- 0.155R 


= 0.219WR + 0.929WH (35) 








Figure 5.1 Passive Earth Pressure on a 180 Degree Central Angle Arch 
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CHAPTER VI 


ANALYSIS OF THE RESPONSE OF THE ARCH TO THE DYNAMIC LOAD 


The arch rib is idealized as a single-degree-of-freedom systen, 
and the flexural resistance and mass of the arch rib per se are ignored. A 
method of dynamic analysis developed by Dr. N. M. Newmark, References (8) 
and (9), is used to predict the response of the arch rib to the given 
pressure pulse. 

For the arches investigated in this thesis, it is found that the 
response of each is in the plastic, or yielded, part of the assumed elasto- 
plastic resistance diagram for the soil. Had the response been less than 
the deflection peeiieed to establish the passive earth pressure, the dif- 
ferential equation of motion stated in the introductory chapter of this study 


would have been used to find the maximum response. 


6.1 Method of Analysis 
For an initially peaked triangular force pulse acting on a single- 


degree-of-freedom system having an elasto-plastic resistance function, 
Figure 6.1, Dr. Newmark presented in Reference (8) the following approximate 


equation to predict the maximum response of the structure: 


Aa 
P 1- = 
eee Vo + ——SR (36) 
qe 1 1+0.7 = 
Me 


where Pp maximum value of the pulse loading 


max 
ca = yield resistance of the structural element 
H = ductility ratio = maximum deflection/yield deflection 


x,/ Xo 





fp 


T = natural period of vibration of the structure = eX \{— 


ty = duration of the loading 


Another technique developed by Dr. Newmark, Reference (9), can 
be used to define the dynamic behavior of the structure for other pulse 
shapes which may be represented by several initially peaked triangular force 
pulses, Figure 6-1. The relationship between the peak pressure and the 


maximum response of the structure is defined approximately as 


a We Pnax [ “1 2 
, ss > eee oe (37) 
u 2 te 1 2 
where Fo F, = damage pressure levels for initially peaked 
triangular force pulses of duration t,/T and 
t,/T, respectively, acting separately 
Ci» Cp = ratio of the peak pressure of the first and 
and second replacement triangles to the total 
peak pressure 
f,,f, = ce, (p_/a,) and c,(p_/4,), respectively, the 


maximum value of the individual force pulses 


To apply the above described methods of dynamic analysis to the 
Situation at hand it is assumed that: 

(1) The configuration of the arch rib is defined by the radial 
deflection of the haunch at a point located at 6/4. The structure is 
idealized as a single-degree-of-freedom system. 

(2) The blast loading on the arch is considered to be a uniformly 
distributed antisymmetrical loading. The unit value of this loading defines 
the forcing function acting on the single-degree-of-freedom system. 

(3) The shape of the blast loading force pulse consists of two 
Maieialiy peaked triangular force pulses. The duration of the first pulse, 


ty» is defined as the time required for the blast wave to engulf the arch 





te 
rib and rise to its maximum pressure at the leeward abutment, a time equal 
to t. plus te The total duration of both pulses is defined as the effective 
duration of the blast loading, ty. 

(4) The resistance provided to the arch rib by the passive earth 
pressure is considered to be a uniformly distributed antisymumetrical load. 
The unit value of this resistance defines the resistance function of the 
single-degree-of-freedom system. 

(5) The weight of the soil associated with the passive earth 
pressure is uniformly distributed along the arch rib. This unit weight 
defines the unit mass of the single-degree-of-freedom systen. 

(6) The passive earth pressure reaches its maximm value when 
the radial deflection of the arch at a point 6/4 is 0.14R. The resistance 
diagram of this pressure is elasto-plastic and defines the resistance func- 
tion of the idealized single-degree-of-freedom systen. 

(7) The period of vibration of the arch is defined as 
T = an mx 7, , where m results fromthe unit weight found in (5), x, is 
O.1%R, and q, is the unit passive earth pressure resistance found in (4). 

The general nature of assumptions (1) through (7) above is show 
in Figure 6.2. 

Numerical values of the parameters in Equation (36) are presented 
in Table VIII. Using these parameters, the maximum responses of arches 
with various radii subjected to 100 and 200 psi overpressures have been 
determined and are plotted in Figure 6.3. It is assumed that the arch is 
buried in soil with a seismic velocity of 2000 fps. The resistance and 
mass Of the arch rib have not been included, and soil arching effects are 


not considered. 





6.2 Sample Calculation 


A sample calculation of the response of an arch under an assumed 


set of conditions follows. 


Let: C. = 2000 fps W= 1M P= 100 psi 
= 2 ft/ms 
w = 120-pef R = 200 in. 
wy lG.67 ft. 
t= at, H, = 0 


The parameters are: 


Paax = 0°37 P,, = 0-37(100) = 37 psi 
oT le 0.06 Po 0.06(100) = 6 psi 
t . R-R sin(? - 8) , R= R sin? 
1 CU 2C 
s 8 
Jct + 
5 
16.67(1 - sin(O - arctan oe 
7 —$—$—$———— ree 4 16.67 
/13 
= 20 ms 
2/3 
_ 100 1/3 
t, = O-4 (= W 
m 
= oo ms 
T = 0.0076VR = 0.0076 V200 = 0.107 sec. 
= 107 ms 
x = 0.001 R = 0.001(200) = 0.2 in. 


aq, = 0.611 wR = aoe 200 = 8.5 psi 


U = 3000 fps 


3 ft/ms 


@ = 180 degrees 





Th 


The force pulse is shown in Figure 6.2. Equations (46) and (37) 
must be evaluated by trial and error. Reference (9), Chart I, may also be 
used to evaluate Equation (36). 

Assume that np = 12. For the first initially peaked triangular 
force pulse having a duration of 20 ms, the damage pressure level necessary 


for a maximum deflection of x = 12x, is 


x 


af 


P 1 = 
max T Jou. 1 Cu 
“qe YH 1 +O) 
q, mt, 1+0.{(T ty) 
i: 
— 
= Box Veh - 1 + 1+ 0-7(107/20) ~ 8.4 = Fy 


For the second initially peaked triangular force pulse having a 
duration of 400 ms, the damage pressure level necessary for a maximim de- 


flection of x = lex is 
m e 


L 
Pnax 10 = oan 1- oT 
“a 7 oo ah - 1 + TP eTMIowoo) = 1:2 = Fo 


The small increase in the peak value of the second pulse due to 
the extension of the pressure line from t = 20 ms to t = O ms, an increase 


of 555 (6) = 0.3 psi, is ignored. 


The actual peak pressure in the first force pulse is 41 psi; the 
actual peak pressure in the second force pulse is 6 psi; and the actual 
ratio of the total peak force to the resistance is 

Prax q 
Bn ee 
oa 8.5 
Equation (357) may now be evaluated to determine if the assumed 


value of up = 12 is satisfactory. 





(5 


Ae) 


Ss 


where c) = 2 = 0.8% Co = gy = 0.16 
p. b.57 ( Set, Oe ee = ~ 1.0 


Therefore the assumed value of 12 for uw is satisfactory, and the 


maximum deflection is 


x = x 
m 7 e 


12(0.2) 


2.4 inches 


6.4 Discussion of Results 

Figure 6.3 contains the quantitative results of the analyses of 
the response of various sized arches. The results may be qualitatively 
summarized as follows: 

(1) For any given overpressure level, the response of 

the arch increases as the radius of the arch decreases. 

This is partially explained by considering how the load 

on the arch and the resistance given to the arch by the 

earth pressure are defined. The blast pressure is de- 

fined as a function of the peak overpressure only. The 

value of the analogous, uniformly-distributed, anti- 

symmetrical load is not a function of the radius of the 

arch. The earth pressure is defined as a function of 

the radius of the arch and varies directly as the radius 


varies. Consequently, the response of the arch to any 





(2) 


(3) 


(4) 
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given overpressure varies indirectly as the radius 
of the arch varies. 

Other factors which probably influence the re- 
sponse of the arch have not been studied in this thesis. 
One of these is the value of the resistance of the arch 
per se. For arches of smaller radii, this resistance 
may be relatively large compared to the resistance given 
the arch by the earth pressure. 

For equal radius arches, the response of the arch 
to any given overpressure level increases as the central 
angle of the arch decreases. 

The effect of the second triangular force pulse is 
very important since it has a relatively long duration. 
Several analyses of the response of a 120 degree arch to 
a 100 psi overpressure were made, in which the effect of 
this second force pulse was not considered. The response 
of a 200 inch radius arch was found to be approximately 
three inches, compared to a value of sixteen inches when 
the effect of the second pulse is included. 

Increasing the depth of cover on the crown reduces 
the effect of the blast load on the arch. The effect of 
the depth of burial of the arch may be related conveniently 
to the ratio of the depth of burial of the arch abutment to 
one-half of the arch span. When this ratio is one, for 
arches having radii of 200 inches or more and central angles 
greater than 120 degrees, the flexural action caused by an 


overpressure of 100 psi resulting from a 1 MT weapon may 





(7 


be ignored. For a 180 degree arch, this ratio is one 
when there is no cover on the crown, and the deflec- 
tion is approximately two inches. For a 120 degree 
arch, this ratio is approximately 0.9 when there is 
a depth of cover of 0.3R on the crown, and the deflec- 


tion is approximately one inch. 


as ‘ 





Parameter 


180 Degree Arch 


H = 0 
Cc 


0.O001R 


1.06wR* 


nm RE 


0.0076 VR 


TABLE VIII 


120 Degree Arch 


H = 0 
Cc 


O.45pP 


m 


Fr 


0.0087 VR 
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PARAMETERS FOR DYNAMIC ANALYSIS OF SOIL RESISTANCE 


120 Degree Arch 
H. = Q.4R 


O.45P 


m 











(a) Single-Degree-of-Freedom 


System 
q(x) 
Gdek-— 
| 
| 
Xe Xm : 
(c) Elasto-Plastic Bilinear 
Resistance 


1 


p(t) 


Pmax 


f, 


(b) Initially Peaked 
Triangular Force Pulse 


p(t) 


Pmax. 


t 
7 te 


(d) Two Initially Peaked 
Triangular Force Pulses 


Figure 6.1 Idealized Single-Degree-of-Freedom System, Force Pulses and 
Elasto-Plastic Bilinear Resistance 





Blast Loading Soil Resistance 


q(x) 


0.00! R Xn 





Soil Resistance Diagram 


Figure 6.2 Idealized Loading and Resistance Functions of an Underground 
Arch Rib Subjected to Blast Pressure 
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CHAPTER VII 


QUALITATIVE DESIGN METHOD FOR ARCH RIBS 


The problem of the structural design of the arch rib has not been 
studied. However, the following qualitative design procedure, based on the 
methods used in this thesis, is suggested: 

(1) Idealize the arch as a single-degree-of-freedom system, the 
configuration of which is defined by the radial deflection of the arch at 
6/4. 

(2) Define the fundamental period of vibration of the structure 
by the mass and resistance of the soil, and the deflection of the arch at 
6/4 necessary to establish the full passive earth pressure. 

(3) Define the forcing function as consisting of two initially- 
peaked triangular force pulses, the first pulse having a duration of te + tw 


the second having a duration of t These forcing functions are uniformly- 


i 
distributed, antisymmetrical, radial loads analogous to the actual blast 
loading. 

(4) Assume a ductility factor uy. The end objective is to 
design the arch rib, not to determine the response of the soil surrounding 
the arch. However, since the arch cannot deflect without the soil ac- 
companying it, this ductility factor also defines the response of the soil. 
The assumed value of uw should be larger than that desired for the actual 
arch rib. 

(5) Evaluate me damage pressure level equation to determine 
the total resistance required of the combined arch rib-soil systen. 

(6) Subtract the soil resistance from the total required 
resistance and design the arch rib with a resistance equal to the remainder. 


The usually recommended allowable stresses for dynamic design are to be used. 
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(7) Determine the yield deflection at 6/4 of the arch section 
designed in (6). The yield deflection mst be stated in the same numerical 
terms as those used in (2) for finding the period of vibration of the 
structure. 

(8) Combine the resistance of the soil and arch rib into an 
elasto-plastic resistance diagram. This is accomplished by making the area 
of the idealized bilinear resistance diagram equal to the area of the 
actual trilinear resistance diagram when both diagrams extend to a numerical- 
ly specified deflection. The ratio x b/ 5. is assumed; the initial slope of 
the elastic part of the actual diagram is maintained; the value of qe is 
adjusted to make the areas equal. Determine a new period of vibration for 
the combined arch-rib soil system, including the resistance and mass of the 
arch section designed in (6). 

(9) Using the elasto-plastic resistance function found in (8), 
evaluate the damage pressure level equation to determine the maximum 
response of the system. The resulting ratio of the maximm deflection to 
the yield deflection of the system is found and compared with the ratio 
x /x, vsed in combining the soil and arch resistances. If these ratios 
are equal, the arch section designed in (6) is satisfactory. If they are 
not equal, the arch section mst be redesigned and (7) through (9) repeated. 

The anticipated results of the design procedure outlined above 
are discussed below without proof. 

There are two principle modes of response of an arch to blast 
loading, Reference (10). One of these is symmetrical, corresponding to 
the compression in the arch, and the other is antisynmmetrical, corresponding 
to the flexural deformation of the arch. Consider an underground arch and 


assume that the arch rib has been designed for the static load of the soil 





cover and the compression load of the blast force. For the arch sizes, 
overpressures, and weapon yield studied in this thesis, it is reasonable 
to expect that such arches would have a yield deflection on the haunch of 
at least one-half inch. For 100 psi overpressure, the maximum deflection 
expected of the haunch, considering the resistance of the soil only, is 
two inches for 180 degree arches with zero depth of cover on the crown, and 
approximately one inch for 120 degree arches with 0.5R cover on the crown. 
Therefore the ratio of the maximum deflection on the haunch to the allovw- 
able elastic deflection of the arch rib does not exceed 4:1, even though 
the reduction of the maximum deflection due to energy absorbed by elastic 
or plastic deformation of the arch rib is disregarded. 

For moderately ductile structures, the above ratio, ductility 
factor u, ranges from 10 to 30, Reference (8). For designing flexural 
members, including reinforced concrete members, the ductility factor usually 
is assumed to be 20. 

It is therefore concluded that if the arch rib is designed so that 
it has a ductility factor of 20 for flexural action, and is covered as de- 
scribed above, no increase in the required resistance for static and blast 
compression load is necessary to accommodate the antisymmetrical mode of 


response for an overpressure of 100 psi. 
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CHAPTER VIII 


CONCLUSIONS AND RECOMMENDATIONS 


8.1 Conclusions 

The following conclusions are drawn from the results obtained 
and the discussions thereof in the preceding chapters. 

(1) Harmonic analysis provides a convenient means of analyzing 
unsymmetrical loads acting on arches. 

(2) The peak values of dynamic loads and their load-time 
curves may be determined by harmonically analyzing the load at suitable 
time intervals. 

(3) The effects of the unsymmetrical blast load actid on the 
arch may be studied conveniently through the use of analogous loads. When 
studying the effect of the blast load on the haunch, the analogous load 
employed is a uniformly distributed radial load which acts in compression 
on the windward side of the arch and in tension on the leeward side. The 
peak value of this load is approximately forty percent of the peak value 
of the overpressure. When studying the effect of the blast load on the 
crown, the analogous load employed is a uniformly distributed radial load 
which acts in compression on the central third of the arch and in tension 
on the outer thirds. The peak value of this load is approximately forty 
percent of the maximum overpressure. 

(4) The effect of the unsymmetrical passive earth pressure load 
on the arch also may be studied through the use of an analogous load. This 
analogous load is a uniformly distributed radial load acting in compression 
on the leeward side of the arch and in tension on the windward side. It 
directly opposes the blast pressure load. Its value varies directly with 


the radius of the arch. 





(5) The peak value of the analogous load employed in studying 
the blast effects on the haunch of the arch increases as the seismic 
velocity of the soil increases. The peak value of the analogous load used 
in studying the crown of the arch decreases as the seismic velocity of the 
soil increases. 

(6) The relative severity of the dynamic load varies inversely 
with the ratio of the horizontal to vertical pressure in the soil. It is 
to be noted, however, that as this ratio increases, there may be a tendency 
for the crown of the arch to move outward. 

(7) The relative rise time of the blast pressure to its peak 
value has little effect on the peak value of the analogous loads. 

(8) The passive earth pressure contributes materially to the re- 
Sistance of the arch. An overpressure of 100 pounds per square inch from a 
one megaton weapon causes a maximum response of approximately two inches 
in the 180 degree arches studied. Subjected to the same overpressure, the 
response of a 120 degree arch with a depth of cover on the crown equal to 
three-tenths of the radius is approximately one inch. This response is 
considered negligible. 

(9) The passive earth pressure developed by the arch decreases as 
the central angle of the arch decreases. It is doubtful if arches with 
central angles of much less than 120 degrees would develop significant 
passive earth pressure under blast loading. 

(10) The rigidity of an arch interferes with the deflection of 
the arch required to develop the passive earth pressure. This suggests that 
flexible arches have certain advantages in blast resistant construction. 

(11) The method of approach presented in this thesis may be 


extended to consider soils with different angles of internal friction, 
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different unit weights, and varying water table levels, as well as to 
cohesive soils, provided the cohesion of the soil can be stated in terms 


of an equivalent angle of internal friction. 


8.2 Recommendations 

It is recommended that the studies made in this thesis be 
extended to: 

(1) Consider how the blast loading cn the arch is affected by 
the soil on the windward side of the arch. 

(2) Consider how the passive earth pressure loading on the arch 
is affected by the increased shear strength of the soil due to the vertical 
blast pressure acting on the leeward side of the arch. 

(3) Compare the arch behavior predicted in this thesis with 
available test data. 

(4) Develop a quantitative procedure for design of the arch rib 
under impulsive loading which includes the effect of the passive earth 


‘ pressure. 
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APPENDIX B 


NUMERICAL CALCULATION RESULTS 


This Appendix contains the results of the numerical harmonic 
analyses made. 

In Tables Bl through BlO are tabulated the vertical pressure-time 
values of blast pressure on a 180 degree arch with zero depth of cover on 
the crown for various combinations of the variables stated in Chapter II. 
Since the 120 degree arch has the same radius, this table also provides the 
vertical pressure-time values for the blast loading on a 120 degree arch by 
using the values between 30 and 150 degrees. 

Tables Bll and Bl2 contain the values of the coefficients p(t) 
for the harmonic components of loading, as found by numerical analysis. 
These values were found by using the vertical pressures in Tables Bl through 
BlO, applying the coefficients from Table I to convert the vertical pressure 
to radial pressure, and then harmonically analyzing these radial pressures. 
These tables also contain the values of the analogous loads p,(t) and p(t). 

The following sign convention is applicable to these loads: 

(a) p(t) is positive when the load acts inwardly 

(b) p(t) is positive when the uniformly distributed load 

aéts in compression on the windward half of the arch 
and in tension on the leeward side 

(c) p(t) is positive when the uniformly distributed load 

acts in compression on the central third of the arch 
and in tension on the outer thirds. 

Originally it was intended to study only twenty-five foot radius 
arches, some with a rise time equal to one-half the transit time and the 


remainder with a rise time equal to the transit time. However, after numerous 


calculations had been made, an error amounting to a miltiple of two was found. 
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In order to use the previous work it became necessary to decrease the arch 
radius to 12.5 feet, and to increase the rise times of some of the arches 
to equal the transit time, and of the remainder to twice the transit time. 
This accounts for the data showing ie = et, - After the error was dis- 
covered, only arches with a twenty-five foot radius and t_ = at, were 
studied. The results are of such a nature that general conclusions may be 
drawn without regard to the disparity in the arch sizes. 

Figures Bl through BO are typical graphs of the values of the 


coefficient p(t) and the analogous loads p,(t) and p(t) - 
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VITA 


Lieutenant Commander Worthen Allen Walls, Civil Engineer Corps, 
U. S. Navy, was born near Little Rock, Arkansas, February 17, 1922. 
Following his graduation from Stephens High School, Stephens, Arkansas, in 
1940, he attended Magnolia A and M College, Magnolia, Arkansas, for two 
years. : 

He entered the Naval Air Corps, U. S. Navy, in June of 1942 and 
after flight training was designated a Naval Aviator and commissioned as an 
Ensign, U. S. Naval Reserve, in 1943. Wartime assignments were primarily 
with the Pacific Fleet and included duty on the West Coast, the Aleutians, 
Hawaii, Saipan, Peleliu, Okinawa, Iwo Jima, and then Japan in the early 
days of the occupation. 

After being released to inactive duty in 1946, he entered the 
University of Arkansas, Fayetteville, Arkansas, from which he graduated, 
with honors, in 1948 with the degree of Bachelor of Science in Civil 
Engineering. He was then employed as an engineer by Phillips Petroleum 
Company, Okmulgee, Oklahoma, until April, 1949. 

In April, 1949, he returned to active duty and was commissioned 
as a Lieutenant (junior grade), Civil Engineer Corps, U. S. Navy. His first 
tour of duty was as Assistant Public Works Officer and Assistant Officer in 
Charge of Construction, Bureau of Yards and Docks Contracts, U. S. Fleet 
Activities, Yokosula, Japan. This duty was followed in 1951 by duty in 
the Public Works Department, U. 5S. Naval Air Station, Corpus Christi, Texas, 
until 1953. 

In May, 1953, he entered the University of Illinois, for graduate 


work in Civil Engineering under a program of postgraduate education 
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sponsored by the U. S. Navy and received the degree of Master of Science 
in Civil Engineering in 1954. 

Thereafter he was assigned as Assistant Resident Officer in 
Charge of Construction, Bureau of Yards and Docks Contracts, Seville, 
Spain, for three years in connection with the Spanish Bases Program. From 
1957 to 1958 he served as Public Works Officer and Officer in Charge of 
Construction, U. S. Naval Ammunition Depot, Concord, California. He 
returned to the University of Illinois in 1958 to continue his advanced 
study in Civil Engineering under Navy sponsorship. 

He is a member of the American Society of Civil Engineers, 


Society of American Military Engineers, and Tau Beta Pi. 
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